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ABSTRACT 

We  propose  mathematical  models  to  describe  the  behaviour  of  smart  obstacles.  In  the  context  of  acoustic 
scattering  a  smart  obstacle  in  an  obstacle  that  when  hit  b  y  an  incoming  acoustic  wave  reacts  circulating 
on  its  boundary  a  “pressure  current  ”  to  pursue  a  given  goal. 

A  pressure  current  is  a  quantity  whose  physical  dimension  is  pressure  divided  by  time.  The  goals 
considered  are:  1)  to  be  undetectable,  2)  to  appear  with  a  shape  and  an  acoustic  boundary  impedance 
different  from  its  actual  ones,  3)  to  appear  in  a  location  in  space  different  from  its  actual  one  eventually 
with  a  shape  and  an  acoustic  boundary’  impedance  different  from  its  actual  ones.  He  mathematical  models 
proposed  for  the  smart  obstacles  are  optimal  control  problem  for  the  wave  equation.  These  optimal 
control  problems  are  studied  analytically  and  solved  quantitatively  using  ad  hoc  numerical  methods. 


1.0  INTRODUCTION 

In  this  paper  a  mathematical  model  for  an  acoustic  time  dependent  scattering  problem  involving  smart 
obstacles  is  formulated.  Smart  obstacles  are  obstacles  that  when  hit  by  an  incoming  acoustic  field  react  in 
order  to  pursue  an  assigned  goal.  The  goal  pursued  by  the  smart  obstacle  considered  in  this  paper  is:  to 
appear  in  a  location  in  space  different  from  its  actual  location  eventually  with  a  shape  and  boundary 
impedance  different  from  its  actual  ones.  We  call  this  goal:  to  appear  as  a  ghost  obstacle.  The  smart 
obstacle  pursues  its  goal  circulating  a  pressure  current  (i.e.  a  quantity  whose  physical  dimension  is 
pressure  divided  by  time)  on  its  boundary.  We  show  that  the  pressure  current  necessary  to  pursue  the  goal 
can  be  determined  as  the  solution  of  a  suitable  optimal  control  problem  for  the  wave  equation. 

The  author  and  its  coworkers  have  studied  similar  models  for  several  other  classes  of  smart  obstacles  in 
acoustic  and  electromagnetic  scattering  (see  for  example  (l)-(6)  and  the  website: 
http://www.econ.univpm.it/recchioni).  The  obstacles  considered  pursue  one  of  the  following  goals: 

1.  to  be  undetectable  (i.e.:  furtivity  problem), 

2.  to  appear  with  a  shape  and  a  boundary  impedance  different  from  its  actual  shape  and  impedance  (i.e.: 
masking  problem), 

3.  to  appear  in  a  location  in  space  different  from  its  actual  location  eventually  with  a  shape  and  boundary 
impedance  different  from  its  actual  ones  (i.e.:  ghost  obstacle  problem). 
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ORGANIZATION 


The  scattering  problems  corresponding  to  1.-3.  have  been  formulated  as  optimal  control  problems  for  the 
wave  equation  (acoustic  case)  or  for  the  Maxwell  equations  (electromagnetic  case)  and  the  first  order 
optimality  conditions  for  these  control  problems  have  been  derived  applying  the  Pontryagin  maximum 
principle  and  solved  with  appropriate  numerical  methods  on  several  test  problems.  The  choice  of  limiting 
the  exposition  to  the  ghost  obstacle  problem  is  motivated  by  the  following  reasons  =  necessity  to  choose  a 
problem  to  fix  the  ideas  and  brevity  and  the  fact  that  the  ghost  obstacle  problem  is  relevant  in  applications 
and  is  harder  than  the  furtivity  and  the  masking  problems.  Several  other  approaches  to  study  smart 
obstacles  have  been  considered  in  the  literature,  see  for  example  (7)-(10). 

2.0  THE  GHOST  OBSTACLE  OPTIMAL  CONTROL  PROBLEM 

Let  QcR^flgcR1  be  two  bounded  simply  connected  open  sets  with  locally  Lipschitz  boundaries 
<9Q,  SQC  and  let  Q  and  Q  g  be  their  closures  respectively.  Let  us  denote  with 
»(*)  =  (wi  Gd  n2  Gd  ni  (*)Y  e  R3 ,  T  e  dQ.  the  outward  unit  normal  vector  to  dfl  in  je  of) .  Since  Q  has  a 
locally  Lipschitz  boundary,  «(x),  x  e  dfl ,  exists  almost  everywhere,  similar  statements  hold  for  the 
outward  unit  normal  vector  to  dflG.  Furthermore  let  QG  be  such  that  Q G  ^0  and  QnQe  =0.  We 
assume  that  Q  and  QG  are  characterized  by  constant  acoustic  boundary  impedances  %  >  0  and  Xo  —  0, 
respectively.  The  case  %  = +co  and/or  Xg  = +00  (he.:  the  case  of  acoustically  hard  obstacles)  can  be 
treated  with  simple  modifications  of  the  formulae  presented  here.  We  refer  to  (Q;^)  as  the  obstacle  and 
to  (D.g;Xg)  as  the  ghost  obstacle.  We  consider  an  acoustic  incident  wave  u‘(x,t\  (x,t) e  R3  x  R, 
propagating  in  a  homogeneous  isotropic  medium  in  equilibrium  at  rest  with  no  source  terms  present  that 
satisfies  the  wave  equation  with  wave  propagation  velocity  c  >  0  in  R3  x  R. 

Finally  we  denote  with  us  {x,t\  (x,/)  e  (r3  \  q)x  R  and  with  uG(x,t),  (x,/)e  (r3  \Qg)xR  ,  the 
waves  scattered  respectively  by  the  obstacle  (Q;^)  and  by  the  ghost  obstacle  (O G ;  x(, )  when  hit  by 
u '  (x,  t),  (x,  t)  e  R3  x  R. 

The  scattered  acoustic  field  us  (x,/),  (x,  t)  e  (r3  \n)xR  is  defined  as  the  solution  of  the  following 
exterior  problem  for  the  wave  equation: 


A  us  (x,t) 


1  d' 


8t 


=  (xq)e (r3  \ fi) 


x  R, 


with  the  boundary  condition: 

-  -^r  (aT)  +  cz  tjt  (x,t)  =  g(x,  t),  (x,  t)  e  do.  x  R, 
at  vn(x) 

where  g(x,t )  is  given  by: 


(1) 


(2) 


g(x,t)  =  {x,t)~  cz  ~~~rx (x,  t),  (x,t)edQx  R, 

at  on\x_) 


(3) 


the  boundary  condition  at  infinity: 


u 


s(x,t)  =  o{^- 


1  /•  — >  +CO,  /  6  R, 


(4) 
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and  the  radiation  condition: 


dus  i  \  1  dus  /  \  ( l) 

- (x, /Jh - [x, t )  =  o\  —  ,  r  — >  +oo, (eR, 

dr  c  dt  \r ) 


(5) 


where  r  =  llxll,  x  e  i?3,  A  =  — -  is  the  Laplace  operator,  c  >  0  is  the  wave  propagation  velocity  and 

' ~l8x f 

0(-)  and  o(-)  are  the  Landau  symbols.  We  note  that  g(x,t),  (xj)edQxR  is  defined  almost  everywhere  and 
that  the  boundary  condition  (2)  can  be  adapted  to  deal  with  the  limit  case  of  the  acoustically  hard 
obstacles,  i.e.  y  =  +oo .  The  obstacle  {0',y)  that  scatters  the  field  us  solution  of  (1),  (2),  (3),  (4),  (5)  is 

called  passive  obstacle.  The  field  usG{x,t),  (xt)g  (r3  \Qg)x  R  scattered  by  the  (passive)  ghost  obstacle 
is  defined  as  the  solution  of  (1),  (2),  (3),  (4),  (5)  when  in  the  problem  defined  above  we  replace  Q.  with 
Qg  and  i  with  yG  ■  Note  that  we  always  consider  the  ghost  obstacle  as  a  passive  obstacle. 

We  consider  the  following  problem: 

Ghost  Obstacle  Problem :  Given  an  incoming  acoustic  field  u‘(x,t),  (x,t)  e  R3  x  R,  an  obstacle  (O;  y)  ■>  a 
ghost  obstacle  (O  G ;  yG )  choose  a  pressure  current  circulating  on  5Q  for  ieR  in  such  a  way  that  the 
wave  scattered  by  [Q',x)  when  hit  by  the  incoming  acoustic  field  u  appears,  outside  a  given  set 
containing  Q  and  QG ,  “as  similar  as  possible”  to  the  wave  scattered  in  the  same  circumstances  by  the 
ghost  obstacle  (fi  G ;  yG ) . 

Remember  that  a  pressure  current  is  a  quantity  whose  physical  dimension  is:  pressure  divided  by  time. 

Our  goal  is  to  model  the  ghost  obstacle  problem  as  an  optimal  control  problem  introducing  a  control 
variable  y/{x,t),  (x,  /  )  e  dQ.  x  R,  that  is  a  pressure  current  acting  on  the  boundary  of  the  obstacle.  To  this 
aim,  we  replace  the  boundary  condition  (2)  with  the  following  boundary  condition: 


(x,  t)  +  cy  (x,  t)  =  g(x,t)  +  (i  +  x)  t),{x,t)edO.Y.  R. 

ot  onyx) 


(6) 


Let  Q  be  a  bounded  simply  connected  open  set  containing  Q  and  QG  with  Lipschitz  boundary  SQ,  and 
let  dsn  ,dsdn  be  the  surface  measures  on  80. s  and  80  respectively. 

We  choose  the  following  cost  functional: 


fA.,m,£(v/)=  \dt  I  j  ^  + x)4us(x,t)-usG(x,t)\  ds8C1<_  +  ^  +  /)iu<;ii/2(x,t)dsdo\. 


en. 


(7) 


an 


where  A,  >  0,  //  >  0  are  adimensional  constants  such  that  A  +  JU  =  1,  and  q  is  a  nonzero  positive 
dimensional  constant.  We  model  the  ghost  obstacle  problem  via  the  following  optimal  control  problem: 

miVcFx++(^)’ 

subject  to  the  constraints  (1),  (4),  (5)  and  (6). 
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This  is  a  legitimate  mathematical  model  of  the  ghost  obstacle  problem.  In  fact  the  minimization  FA  e 

makes  small  us -usG  for  (x, i) e dO.E  x R ,  that  is  makes  small  us -usG  for  (x,?)e(r3 \Qe )xR  and  makes 

small  the  “size”  of  the  pressure  current  used  while  the  constraints  (1),  (4),  (5),  (6)  guarantee  the 
satisfaction  of  the  dynamic  conditions  associated  to  the  problem  considered. 

The  set  C  is  the  space  of  the  admissible  controls  that  we  leave  undetermined.  The  obstacle  (Q;  %)  that 

generates  the  scattered  field  u  solution  of  (8),  (1),  (4),  (5),  (6)  is  called  smart  or  active  obstacle. 

Note  that  in  (7)  the  choice  QG  c  Q,  Q,  =  Q  gives  the  masking  problem  and  that  the  choice  flG  =  </), 

Q£  =Q  gives  the  furtivity  problem. 

3.0  THE  FIRST  ORDER  OPTIMALITY  CONDITIONS 

Let  us  make  the  following  assumptions:  let  be  the  usual  spherical  coordinate  system  in  R3  with 

pole  in  the  origin,  let  B  be  the  sphere  with  center  the  origin  and  radius  one  and  let  dB  be  its  boundary,  we 
assume  that: 

(a)  the  boundary  of  the  obstacle  Q  is  a  starlike  surface  with  respect  to  the  origin,  that  is  Q  and  8B  can  be 
represented  as  follows: 


Q  =  {x  =  rx  e  R3  0  <  r  <  e  dB], 

(9) 

dd  =  {x  =  rx  e  RJ  \r  =  <^(x),x  e  8B }, 

(10) 

where  <%(x)>  0,x  e  dB ,  is  a  single  valued  function  defined  on  dB  that  is 

regular  for  the  manipulations  that  follow; 

the  sets  Q,  and  GQt.  can  be  represented  as  follows: 

assumed  sufficiently 

=  {x  =  rx  e  R3|0  <  r  <  (^(x)+s),x  e  8b\s  >  0  > 

(11) 

8Qe  =  {x  =  rx  e  R3|r  =  £(x)+  s,x  e  dB}, e  >  0  ■ 

(12) 

for  a  suitable  choice  of  s  >  0. 

Note  that  the  assumptions  (a)  and  ( b )  are  only  one  of  many  other  possible  choices  of  assumptions  that  can 
be  made  to  guarantee  the  satisfactory  solution  of  the  model  (8),  (1),  (4),  (5),  (6).  This  choice  is  made  just 
to  fix  the  ideas  and  to  keep  the  exposition  simple. 

Under  the  assumptions  (a)  and  ( b ),  applying  the  Pontryagin  maximum  principle  the  optimal  state 
trajectory  us  and  the  corresponding  adjoint  variable  trajectory  (p  satisfy  the  necessary  first  order 
optimality  conditions  associated  to  the  optimal  control  problem  (8),  (1),  (4),  (5),  (6),  that  is  they  are  the 
solution  of  the  following  exterior  problem  for  a  system  of  two  coupled  wave  equations: 


Aus(x,t)—  \  d  u ^  (x,t)  =  0,  (x,?)e  (r3  \q)xR, 


c  dt 


( x,t)  =  0 


r  — >  +oo,  teR 


\rj 


(13) 

(14) 
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dus  /  \  1 dus  /  \  fl) 

- \x,t)+ - [x,t)  =  o\  —  L  r— >+  co,  ?eR> 

dr  c  dt  l  rl 


dus  /  \  dus  /  \  [  \ 

--z-{x,t)  +  cz—n-(x,t)=g(x,t)- 
dr  dn\x) 


il  +  x) 

2  ns 


<p(x,t),(x,t) e  ofi x  R  ■ 


A  <p(x,t)~  \  ^  f  (x,t)  =  0,  (x,t)  e  (r3  \q)xR, 
c  dt ' 


^>(x,t)=  (9 


rr 

UJ 


r  — »  +oo,  JgR. 


—  (x,r)+— — (x,t)  =  of— \  r  — >  +oo,  (gRi 
Sr  '  c  St  '  UJ 

’  dT (-’*)- CZd«(x)(-)=  -241+x)//g]Wx+^,^-<fx+£Q>?|LMee£J><R, 


lim  us  (x,  t)  =  0,  x  e  R 3  \  fi  . 


lim<^(x,t)  =  0,xeR3  \  Q, 

t— >+00 


(15) 

(16) 

(17) 

(18) 

(19) 

(20) 

(21) 

(22) 


where  fE (x/  x  Lx  e  8Q  is  the  function  defined  by: 


fE  I  |j=j-  I  =  fe  (x(&,  $))  =  6 ’f) ,  x  g  8Q.,  x  =  2L  g  SB, 


u{U, 


(23) 


0  <0<  7T,  Q<(j)  <2ji  , 


v(0,0)  =  fj |  ff  |  sin2  9  + 

,66 


0<6<7r,0<(j)<27r, 


dt; 


+  £2  sin2  9, 


(24) 


v*(0>0)  =  (f  +  £)](^  1  sin26l  +  [^f]  +(4  +  e)2  sin2  6, 

0  <  9  <  ;r,0  <  </>  <  2k  ■ 


(25) 


The  relation  between  cp  and  the  optimal  control  y/  solution  of  problem  (8),  (1),  (4),  (5),  (6)  is  the 
following  one: 


y/{x,t)  = - —  <p{x,t\{x,t)e  3QxR. 

2  pq 


(26) 
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Let  us  point  out  that  we  have: 

dssn  =  v(0,(f)d6d(j),0  <6<  <  (j)  <  2n , 


(27) 


and 


dsme  =  us  {0,  (f)dOd(j)fi  <0<  n,0  <  (f>  <  2n , 


(28) 


4.0  NUMERICAL  SOLUTION  OL  THE  EXTERIOR  PROBLEM  (13)  -  (22) 

Numerical  methods  to  solve  the  exterior  problem  (13)-(22)  have  been  developed  in  [5],  [6].  These 
methods  belong  to  the  class  of  the  operator  expansion  methods  and  are  highly  parallelizable.  Some 
numerical  experiments  proving  the  validity  of  the  control  problem  proposed  as  mathematical  model  of  the 
ghost  obstacle  problem  are  shown  in  the  website  http://www.econ.univnm.it/recchioni/wl  1 . 


5.0  EXTENSION  AND  CONCLUSIONS 

The  work  presented  can  be  extended  to  a  new  class  of  smart  obstacles  that  pursue  the  following  goal: 

4.  one  of  the  goals  specified  in  the  Introduction  restricted  to  a  definite  band  in  the  frequency  space. 

For  reasons  of  brevity  we  restrict  our  attention  to  the  definite  band  ghost  obstacle  problem  in  the  acoustic 
case.  This  problem  is  formulated  as  an  optimal  control  problem  for  the  wave  equation.  The  acoustic 
definite  band  masking  problem  and  furtivity  problem  can  be  treated  similarly.  We  consider  the  study  of 
these  problems  as  preliminary  to  the  study  of  the  corresponding  problems  in  the  electromagnetic  case 
where  the  wave  equation  must  be  replaced  with  the  Maxwell  equations. 

We  note  that  restricting  the  goal  pursued  to  a  definite  band  in  the  frequency  space  modifies  substantially 
the  mathematical  formulation  of  the  problems  under  scrutiny.  In  fact  the  optimal  control  problems  used  to 
model  problems  1),  2),  3)  in  particular  the  cost  functionals  that  must  be  minimized  in  order  to  model 
appropriately  the  problem  formulated  only  on  the  desired  band  in  the  frequency  space  are  “nonlocal”.  That 
is  the  presence  of  the  definite  band  makes  necessary  the  use  of  suitable  convolutions  involving  the  anti 
Fourier  transform  of  the  characteristic  function  of  the  definite  frequency  band  in  the  definition  of  the  cost 
functional.  Consequently  the  first  order  optimality  conditions  of  these  new  optimal  control  problems 
change  substantially  and  cannot  be  deduced  from  those  derived  in  [2],  [3],  [5]  and  here  for  the  optimal 
control  problems  1),  2),  3).  That  is  the  first  order  optimally  conditions  are  not  expressed  by  two  wave 
equations  coupled  by  local  boundary  conditions  as  in  [2],  [3],  [5]  and  here  but  the  coupling  between  the 
two  wave  equations  is  given  by  nonlocal  (in  time)  boundary  conditions.  As  a  consequence  the  way  of 
solving  the  first  order  optimality  conditions  must  be  changed. 

We  can  conclude  that  the  idea  of  modelling  the  smart  obstacles  using  optimal  control  problems  is  an 
interesting  idea.  Moreover  the  work  developed  until  now  with  the  model  proposed  can  be  profitably 
extended  in  several  directions  such  as  the  study  of  closed  loop  controls,  finite  horizon  controls,  or  the 
study  of  inverse  problems  involving  smart  obstacles.  These  are  challenging  mathematical  questions  whose 
solution  can  be  very  valuable  in  practical  applications. 
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Smart  (or  active)  obstacles  are  obstacles  that  when 
illuminated  by  an  incoming  field  react  actuating  a  policy 
in  order  to  pursue  an  assigned  goal.  The  design  of  smart 
obstacles  can  be  improved  by  the  availability  of 
satisfactory  mathematical  models  of  them.  We  propose 
the  use  of  models  based  on  optimal  control  problems  to 
describe  the  behaviour  of  smart  obstacles. 

Let  us  restrict  our  attention  to  the  context  of  time 
dependent  and  time  harmonic  acoustic  or 
electromagnetic  scattering.  In  this  context  the  smart 
obstacle  in  order  to  pursue  its  goal  circulates  on  its 
boundary  a  pressure  current  (acoustic  case)  or  a  surface 
electric  current  density  (electromagnetic  case). 
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The  goal  pursued  by  the  smart  obstacle 
considered  in  our  work  is  one  of  the  following: 

to  be  undetectable  (furtivity  problem), 
to  appear  with  a  shape  and  boundary  impedance 
different  from  its  actual  ones  (masking  problem), 
to  appear  in  a  location  in  space  different  from  its 
actual  one  eventually  with  a  shape  and  boundary 
impedance  different  from  its  actual  ones  (ghost 
obstacle  problem), 

one  of  the  previous  goals  restricted  to  a  given  band  in 
the  frequency  space  (definite  band  problems). 

These  problems  are  translated  in  optimal  control 
problems  for  the  wave  equation  (acoustic  case)  or  for 
the  Maxwell  equations  (electromagnetic  case). 
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We  consider: 


1) Direct  scattering  problems 

In  this  case  the  optimal  control  problems  associated  to  the 
smart  obstacles  give  a  way  of  characterizing  and 
computing  the  currents  needed  to  pursue  the  assigned  goal 
as  optimal  solutions  of  the  mathematical  problems 
considered, 

2)  Inverse  scattering  problems 

That  is  starting  from  some  knowledge  of  Mfar  field  data” 
generated  by  the  smart  obstacle  in  correspondence  of 
several  known  incoming  fields  and  the  knowledge  of  the 
’’nature"  (furtive,  masked  ghost)  of  the  smart  obstacle  we 
want  to  reconstruct  the  shape  of  the  obstacle.  In  this  case 
the  optimal  control  model  gives  the  mathematical 
framework  where  to  formulate  the  inverse  problem. 

The  work  presented  is  joint  work  with: 

L.  Fatone,  F.  Mariani,  G.  Pacelli,  M.C.  Recchioni. 
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Direct  Scattering  Problems 

See  the  papers: 

F.  Mariani,  M.C.  Recchioni,  F.  Zirilli:  "The  use  of  the  Pontryagin 
maximum  principle  in  a  furtivity  problem  in  time  dependent 
acoustic  obstacle  scattering",  Waves  in  Random  Media  JJ_,  (2001), 
549-575. 

L.  Fatone,  M.C.  Recchioni,  F.  Zirilli:  "A  masking  problem  in  time 
dependent  acoustic  obstacle  scattering",  ARLO-Acoustic  Research 
Letters  Online  5,  Issue  2,  (2004),  25-30. 

L.  Fatone,  M.C.  Recchioni,  F.  Zirilli:  "Furtivity  and  masking 
problems  in  time  dependent  electromagnetic  obstacle  scattering", 
Journal  of  Optimization  Theory  and  Applications,  121,  (2004),  223- 
257. 

L.  Fatone,  G.  Pacelli,  M.C.  Recchioni,  F.  Zirilli:  "The  use  of 
optimal  control  methods  to  study  two  classes  of  smart  obstacles  in 
time  dependent  acoustic  scattering",  to  appear  in  Journal  of 
Engineering  Mathematics, 
and  the  websites: 

http://www.econ.univpm.it/recchioni/w6 
http://www.econ.univpm.it/recchioni/wlO 
http://www.econ.univpm.it/recchioni/wl  1 

The  websites  contain  animations  and  virtual  reality  applications 
representing  data  obtained  in  numerical  experiments. 
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Inverse  Scattering  Problems 

See  the  papers: 

L.  Fatone,  M.C.  Recchioni,  A.  Scoccia,  F.  Zirilli: 
“Direct  and  inverse  scattering  problems  involving  smart 
obstacles”  Journal  of  Inverse  and  Ill  Posed  Problems  13, 
(2005),  247-257. 

L.  Fatone,  M.C.  Recchioni,  F.  Zirilli:  “A  method  to 
solve  an  acoustic  inverse  scattering  problem  involving 
smart  obstacles”  submitted  to  Waves  in  Random  and 
Complex  Media. 

and  the  website: 

http://www.econ.univpm.it/recchioni/wl3 

The  website  contains  animations  and  virtual  reality 
applications  representing  data  obtained  in  numerical 
experiments. 


6 


For  a  general  reference  to  our  work  on  scattering 
see: 

http://www.econ.univpm.it/recchioni 
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We  treat  furtivity,  masking,  ghost  obstacle  and 
definite  band  problems  in  acoustic  and 
electromagnetic  obstacle  scattering  in  a  similar  way. 

For  simplicity  in  this  talk  we  begin  concentrating  on 
the  direct  problem  in  the  simplest  case:  furtivity 
problems  in  acoustics. 

We  consider  a  furtivity  problem  in  the  context  of  time 
dependent  three  dimensional  acoustic  obstacle 
scattering. 
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3 

Let  Q  cz  R  be  a  bounded  simply  connected  open 
set  with  locally  Lipschitz  boundary  d€l  and  let  Q 
be  its  closure.  Let  n_(x)  be  the  outward  unit  normal 
vector  to  d  Q  in  x  e  d  Q . 

We  assume  that  there  exists  a  >  0  such  that: 

5<2=|_xeR3  x<ajc:Q. 

We  consider  acoustic  waves  propagating  in  a 
homogeneous,  isotropic  medium  in  equilibrium  at 

'i 

rest  filling  R  \  Q  with  no  source  terms  present. 

The  obstacle  Q  is  the  scatterer  that  we  assume  to  be 
characterized  by  a  constant  boundary  acoustic 
impedance  Z>0. 

The  limit  cases  of  acoustically  soft  obstacles  (x  =  0) 
and  acoustically  hard  obstacles  (z  =  +00)  can  be 
considered. 


Let  u*  (x,t\(x9t)e  R3  x  R  be  the  incoming  acoustic 
wave,  we  assume  that  ul  (x^)  satisfies: 


At/ i (x, ^ ) -  ^  \  (x,t)=  0,  (i,/)gR3xR 

c2  dt2 
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is  the  wave 


3 

where  ^  _  X  ~Z  2  and  c  >  0 

i= 1 

propagation  velocity. 


When  the  incoming  wave  hits  the  obstacle  Q  a 
scattered  wave  us  (x,t\(x,t)  e  (r5  \q)xR  is  generated. 


the 


For  "passive"  obstacles  us(x,t)  is  defined 
solution  of  the  following  problem: 


s  1  d^us 


(1)  Au  -7^ 


0,  (x,/)g(r^\q)xR 


3  u  3  u  /  \  /  \  _ 

(2)  -——  +  cz-T—  =  g(x,t)i(x,t)ednxR9 


s,  \  m 

(3)  (x,r)  =  0  -  ,  r  = 


x  — >  coj  e  R 


✓  w5  1  d  us  ( l'] 

(4)  — —  +  — —  =  0  -  ,r^co,  reR, 

c  V  r 


where 
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g{x,t) 


dul  d  ul 

- CX — 

d  t  d  n 


(x,^)g^QxRi 


The  (direct)  furtivity  problem  studied  here  can  be 
stated  as  follows:  given  the  incoming  wave  packet 

ul ,  the  obstacle  Q  and  its  acoustic  boundary 
impedance  %  choose  a  control  function,  in  a  suitable 
class  of  admissible  controls,  in  order  to  minimize  a 
cost  functional  that  roughly  speaking  measures  the 
"magnitude"  of  the  scattered  wave. 

The  control  function  is  defined  on  the  boundary  of 
the  obstacle  d  Q  for  t  g  R . 


The  presence  of  this  control  function  changes  the 
nature  of  the  obstacle  from  being  a  "passive" 
obstacle  in  being  an  "active"  obstacle,  that  is  a 
"smart"  obstacle. 
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Let  C  be  the  space  of  the  admissible  controls,  we 
want  to  solve  the  following  optimal  control  problem: 


min 

i//<=C 


subject  to: 

(6)  Xus  — \~~Y~  =  0’  Me(Rivn)xR, 

c  d  i 


cus 
d  t 


+  c  X 


dus 

dn 


g(x, t)  +  (l  +  x)  vi*, t\ (i,^)e^QxR? 


us  (x,t)= d 


\r  J 


r  oo,/  e  R 


dus  1  d  us  ( O 
(9)  —  + - =  * 


d  r  c  d  t 


,  r  —>  oo,  t  e  R 


\rj 


where: 

(10 )Fx,ju(v/)=  \dt  \{i  +  x)  [d'l  + 


R  d£l 


ds 


dfl 


A>0,  ju>0,A  +  ju  =  \  and  ^  >0  is  a  dimensional 
constant. 

Note  that  A,  fd  are  adimensional. 
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Remind  that: 


Fx ,//W=  \dt  K1  +  t) 

R  dQ. 


+  dsd  a . 


j 


Remarks: 

s 

1 .  u  depends  on  \j/  via  the  boundary  condition  (7), 

2.  2  =  0,  //  =  1  corresponds  to  the  "passive" 
obstacle,  in  fact  in  this  case  y>  =  0  minimizes 

FX,nW), 

3.  2  =  1,  //  =  0  corresponds  to  the  "undetectable" 

1 

obstacle,  in  fact  in  this  case  ¥  =  ~~ - g  gives 

l  +  z  s 

us  =  0  and  Fx,n (ft) =  0 .  That  is  the  proposed 
choice  of  Y  minimizes  FX,ju  (y)  and  makes  the 
obstacle  undetectable  (i.e.  us  =0). 

The  functional  (y/)  generates  interesting 
problems  when  A>0,  //  >  0, 2  +  //  =  1 . 

In  this  case  minimizing  (10)  we  make  us  "small"  in 
(r^  \  q)x  R  (we  make  the  obstacle  furtive).  The 
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explicit  presence  of  \j/  in  (10)  (i.e.  ju>  0),  implies 
that  a  cost  must  be  paid  to  achieve  furtivity. 

The  assumption  that  the  control  function  \j/  is 
defined  on  ^QxR  is  a  natural  assumption  in  many 
applications. 

Dimensionally  \j /  is  pressure  divided  by  time  so  that 
we  call  \\f  "pressure  current". 

Let  us  apply  the  Pontryagin  maximum  principle  to 
the  previous  optimal  control  problem. 

Note  that  \q)xR  is  the  "state  trajectory" 

and  i//(x,t\(x,t)ed Q  x  R  is  the  "control  function". 

We  begin  choosing  C  the  set  of  the  admissible 
controls  and  the  space  U  that  is  supposed  to  contain 

US  I^Q. 


We  define: 

D+  -  'Ir3  \  q)x  R  — »  R  such  that  F 


dF 

dn 


d  QxR 


dF 
d  t 


d  QxR  G  L2  (d  Q  x  R), 


<^QxR  , 


F  satisfies  the  wave  equation  (1),  the  boundary 
condition  at  infinity  (3)  and 


dF  -  1  d  F 

- 1 - o 

dr  c  d t 


f-1 

UJ 


,  r  — >  oo,  t  <e  R 


C  =  {/  eL2 (i?Qx  R),  f(x, t)<=  Lx(dQ),t  e  R 


such  that  3  F  e  D_  such  that 


^I^QxR  =  /,  lim  F(x,t)=  0,xeRJ\Q 
/■— >-00 


U  =  {f  gL2 (<?  a X  R), /(x,t)  g  L00 (<?  Q),(gR 


such  that  3  F  eD+  such  that 


F  dflxR  = /’  F(x,t)=0,xeRJ  \  C1 

t  — >— GO 


Note  that  C  and  U  are  vector  spaces. 

The  Hamiltonian  associated  to  the  optimal  control 
problem  considered  previously  is  given  by: 

h(us ,y/,<p,<p0,t)=  jds#n(x)  |  <p0(x,t)  (l  +  x) 

dQ.  y 


't  (us(x,t)f  +  jU^(i//(x,t)f  +  <p(x,t )  cx^-{x,t)-{i  +  x)¥{x,t)~  g{x,t) 


d  n 


The  functions  cp(x, t\ cp0  (x^)  are  called  ’’adjoint 
variables". 

Note  that  we  can  choose  (p0  =  . 

Let  w!(x,4^)eRjxR  be  an  incident  wave  such 

that  V^gC  the  corresponding  us  solution  of  (6), 
(7),  (8),  (9),  satisfies  the  condition: 

lim  us(x,t)  =  0,  x  e  \Q 
t— >-00 


d  t 


d 


S  u  u  C/  u  n  \ 

and  u  ^qxr>^t— <?QxR>^^  ^QxR>“  w £2xKf 


dn 


Using  the  Pontryagin  maximum  principle  it  can  be 
shown  that  the  optimal  solution  of  the  optimal 
control  problem  (5),  (6),  (7),  (8),  (9)  satisfies  the 
following  set  of  equations: 


(ii) 


Air  - 


1  d2us 

77T 


(12  )us(x,t)  =  0 


f-1 

UJ 


=  0,  (i,()e(R3\n)xR) 
r  — »  oo  t  e  R 

? 


(13) 

(14) 


d  us 

- 4 

d  r 


1  dus 
c  d  t 


,  r  — »  oo,  t  e  R 


dus  dus 

- 4cr - 

d  t  d  n 


=  g 


^  +  ^  (p,  (x,  t)  g  d  Q  x  R 
2/u%  V  ' 
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2 

(15)  Ap- 4^4  =  0,  Mef^VnjxR, 

c1  dt1 

(16)  = 


r  ^ooJeR 


(17) 

(18) 


\rj 

d  (p  \  d  cp  _  (\ 


dr  c  d  t 


d  (p 


o 


—  |,  r  — »  oo,  t  e  R 
VT 


=  -2A(l+r)  us  ,(x,t)  e  d  Qx R, 


(19)  lim  u 
t—>—  00 


lim  (p  -  0,  x  e  \  f2 . 

>+C0 


The  relation  between  cp  and  the  optimal  control  V  of 
problem  (5),  (6),  (7),  (8),  (9)  is  the  following  one: 

(20)  i/>(x,t)  =  --]—(p(x,t),(x,t)EdQxR  0  <  pi  <  1. 

That  is  the  solution  of  the  optimal  control  problem 
(5),  (6),  (7),  (8),  (9)  can  be  characterized  through 
(20)  via  the  solution  of  an  exterior  problem  for  two 
coupled  wave  equations,  that  is  problem  (11),  (12), 
(13),  (14),  (15),  (16),  (17),  (18),  (19). 

This  fact  makes  possible  to  solve  the  control 
problem  for  the  smart  obstacle  approximately  at  the 
same  computational  cost  than  the  cost  necessary  to 
solve  the  scattering  problem  for  the  passive  obstacle. 
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Let  us  consider  now  masking  problems. 


Acoustic  Masking  Problem 

Given  the  incoming  acoustic  wave  packet  u\  the 
obstacle  Q  and  its  boundary  impedance  X ,  and 

given  an  obstacle  D  such  that  D  cz  Q ,  with 

boundary  impendance  Xa ,  choose  a  pressure  current 
(i.e.  a  control  function)  defined  on  the  boundary  of 
the  obstacle  G,  d  Q ,  for  te  R  in  order  to  minimize  a 
cost  functional  that  measures  the  "difference" 
between  the  wave  scattered  by  Q  and  the  wave 

scattered  by  D  when  hit  by  ul . 


* 

Note  that  the  couple  D,  Xa  is  called  "the  mask". 


Let  usD  be  the  scattered  field  generated  by  D,  Xa 
when  hit  by  u1  (D  cz  Q ). 

The  masking  problem  considered  can  be  treated  as  a 
furtivity  problem  with  the  cost  functional  given  by: 


(21)  F2A*iM=  lRdtf  n(1+x)MuS  -uSd)  2+li^J  . 


Inverse  Scattering  Problems 

We  limit  our  attention  to  time  harmonic  inverse 
scattering  problems. 

7  _  co  s 

Let  K~  —  and  uk,a  be  the  time  harmonic  field 

C  — 

u  S  t ),  (x,  t )  e  (r ^ 3  \  Q )x  R  scattered  by  the  smart 
obstacle  when  hit  by  the  time  armonic  plane  wave 

uk,a  =  similarly  let  <Pk,a  be  the  corresponding 

time  harmonic  component  of  the  adjoint  variable 

^(x^),(x,^)e[r3 \q)xR  and  uD,k,a  be  the 
corresponding  time  harmonic  component  of  the  field 
scattered  by  the  mask  usD  (x,t)  e  (r3  \  q)x  R . 
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Arguing  as  before  it  can  be  seen  that  for 

$ 

k  e  R,<2  e  d  B  the  functions  uk,a'>(Pk,a  satisfy  the 
following  exterior  boundary  value  problem: 

[Aul  a  +  k2uk,a  )(*)  =  x  e  R3  \  Q, 

(An-,«  +  k2(Pk,a  )(*)  =  0,  x  e  R3  \  a. 


(22) 

(23) 


(\-X)ik(pkJx)-{\-X)x 


+ 


=  0,  x  e  dQ, 


with  the  following  conditions  at  infinity: 


(26) 

(27) 
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where 

bk  a (x)  =  ~ikelk(-’-\l  +  z(n(x\  &)\  *  G 

Let  uk,a  (— )’  %;,«(*)>*  G  R3  be  the  solution  of 

(22)-(27).  The  function  %;,«(*)’  * G  ^  ^ , 
solution  of  (22)-(27)  is  an  auxiliary  unknown  related 

to  the  optimal  control  vi*)  =  x  g  5Q  ?  by  the 
following  relation: 


\<Pk,a(x\  xedQ 
€  ~ 


Due  to  the  first  order  optimality  condition  (22)-(27) 
the  solution  of  the  acoustic  time  harmonic  direct 
masking  problem  consists  in  finding  the  scattered 

acoustic  field  uk,a  and  the  auxiliary  function  Pk,a 
solutions  of  (22)-(27)  and  in  determining  the  optimal 
control  function  via  equation  (28). 
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•  s  •  • 

In  particular  uk,a  is  solution  of  the  Helmholtz 

-5 

equation  (22)  in  R  \  Q  and  satisfies  the 
Sommerfeld  radiation  condition  (26)  at  infinity,  this 
implies  that: 


u 


s 

k,a 


ikr 


r 


Fu  (x,  k,a)+0 


m 

\r2  j 


r  —>  oo 


The  term  Fu{x,k,a)  [s  called  "far  field"  of  uk,a  . 

Similarly  k, a  (T  =  complex  conjugate)  is  solution 

of  the  Helmholtz  equation  (23)  in  R  \  H  and 
satisfies  the  Sommerfeld  radiation  condition  (26) 
((26)  complex  conjugate),  this  implies  that: 


<P 


r  —>  oo 


The  term  ^  (x,k,a)  called  "far  field"  of  Vk,a  . 
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Acoustic  Time  Harmonic  Masking  Inverse 
Scattering  Problem 

From  the  knowledge  of  several  far  fields  generated 
by  the  masked  obstacle  when  hit  by  known  incident 
acoustic  waves  when  the  optimal  pressure  current  is 
active  and  of  the  acoustic  boundary  impedances  of 
the  obstacle  and  of  the  mask  find  the  shape  of  the 
obstacle. 


23 


We  generalize  to  the  case  of  smart  (masked) 
obstacles  the  Herglotz  function  method  used  to 
reconstruct  the  shape  of  an  obstacle  from  far  field 
data  in  the  case  of  passive  obstacles. 

This  method  has  been  introduced  in: 

D.  Colton  and  P.  Monk,  “The  numerical  solution  of 
the  three  dimensional  inverse  scattering  problem  for 
time  harmonic  acoustic  waves”,  SIAM  Journal  on 
Scientific  and  Statistical  Computing,  8,  (1987),  193- 
200,  and  developed  by  several  authors  including: 

L.  Misici,  F.  Zirilli,  “Three-dimensional  inverse 
obstacle  scattering  for  time  harmonic  waves:  a 
numerical  metnod”,  SIAM  Journal  on  Scientific 
Computing,  15,  (1994),  1174-1189. 
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Let  E  be  the  set  of  the  eigenvalues  of  the  matrix 

(A  (A 

Laplace  operator  A  =  |^0  a J  ac^ng  on  two  dimensional 

vector  functions  defined  in  the  domain  Q  with  the 
homogeneous  boundary  condition  corresponding  the 
to  boundary  condition  (32)  on  dd  . 


2 

Definition  1  Given  -  k  let  W_k  be  the  unique 

solution  of  the  equation: 


{*Wk)(y)+k2Wk(y)=0,  yeQ 


(31) 


with  the  boundary  condition: 
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(32) 
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where 


We  say  that  Q  is  a  smart  Herglotz  domain  at  wave 
number  k  with  respect  to  the  boundary  condition 

(32)  if  there  exists  f^(*)  =  G  such  that: 


rSk,u  (*)) 


x  e  dB , 


and 

Kk (z)  =  ds8B (x)e'k^’-gk (Z),  yea 


It  is  easy  to  see  that  the  class  of  the  smart  Herglotz 
domains  at  wave  number  k  is  not  empty  in  fact  for 

example,  when  -  k 2  the  passive  sphere  is  a 

smart  Herglotz  domain. 
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Let  us  define: 


'Foy(x,k,af 

Fot7p(x,k,a)/ 


,  x  g  dB,  «  g  dB , 


and  for  a,b  e  R  let 
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Now  we  can  describe  the  numerical  method  that  we 


propose  to  solve  the  inverse  scattering  problem 
involving  smart  obstacles  considered  here. 

Let  Nl,N2,N3,N4,  be  four  positive  integers,  let 

Xi=  K'  >0I ,*  =  1,2,...,^}  be  the  set  of  the 
wave  numbers  considered,  X2  =  \aie  dB,  i  =  1,2,...,  N2 } 
be  the  set  of  propagation  directions  of  the  incoming 
waves  considered,  GdB,  i  =  \,2,-;Nl}  be  the 

set  of  directions  where  the  far  fields  generated  by  the 
incoming  waves  are  measured  and 
Yj4=  & edB,  i  =  1,2,...,n4}  be  a  second  set  of 

directions  where  the  far  fields  generated  by  the 
incoming  waves  are  measured.  Moreover  let 

X3  and  X4  be  such  that  X3nX4=0- 


Let  n; i  be  a  positive  integer  and  let  us  choose 
A, (o,l),  /  =  1,2..., na  such  that  <2/+1,  i  =  l,2,...,NA -1? 
the  numerical  algorithm  that  we  propose  to 
reconstruct  dQ  is  based  on  the  following  iterative 
procedure:  let  i  =  1 

zo )  set  2  =  A- ; 

h )  choose  an  initial  approximation  d  Q  *  of  the 
boundary  dQ; 

z2 )  using  the  current  approximation  dQ*  of  dQ 
from  the  knowledge  of  the  value  given  to  the 
parameter  X  appearing  in  the  first  order 
optimally  conditions  (22)-(27)  and  the 

knowledge  of  the  measures  of  the  far 

fields  ,u  {x,k,a),k  e  Z3  u  Z4 

approximate  the  far  field 

F0  ,  (x,k,a),Fo,p(x,k,a),k  e  Zlsa  e  Z2,xg  dB  ; 

z3 )  from  the  relation: 

aT=  f ds8B (*) (f(v ki , «7- ),rk^b (*))>  *  = 

ki  dB  '  (36) 

j  =  1,2 ,...,N2,a,b  e  R, 

and  from  the  knowledge  of  the  far  fields 
associated  to  and  (pka  reconstruct  an 

approximation  of  the  vector  Herglotz  kernels 
Sk  ’  i  =  associated  to  the  smart  Herglotz 

domain  Q; 


U)  from  equation  (34)  and  the  knowledge  of  the 
approximate  Herglotz  kernels  reconstruct  the 
associated  vector  Herglotz  wave  functions 
Wki,  i  =  1,2,..., Ni; 

h )  from  equation  (32)  and  from  the  knowledge  of 
the  vector  Herglotz  wave  functions 
?  z  =  l,2,...,Ah  reconstruct  a  new 

approximation  d Cl  of  the  boundary  of  the 
smart  obstacle; 

h)  if  the  “distance”  between  the  old  and  the  new 
approximation  of  dQ,  that  is  the  “distance” 

between  dQ*  and  dQ  is  “small”  set 

dat  =dCl  and  go  to  Step  h)  otherwise  set 

d  Q*  =  d  Q  and  go  to  Step  h); 

ll)  if  i  =  N; i  stop  otherwise  set  i  =  i  +  l  and  go  to 
Step  zo). 

The  reconstruction  of  the  couple  (T,Q)  chosen 
among  the  solutions  i  =  l,2,...,NA  obtained 

with  the  previous  procedure  is  the  one  where  the 
boundary  condition  (32)  is  satisfied  more  accurately. 
That  is  if  the  boundary  condition  (32)  is  imposed  in 
the  least  squares  sense  the  one  that  minimizes  the 
remainder  at  the  end  of  the  least  squares  procedure.  30 


The  first  experiment  shown,  involves  a  furtive 
acoustically  soft  (i.e.  x  =  0)  double  cone  (see  Fig.  1). 
The  double  cone  is  obtained  placing  two  cones  of 
the  same  height  and  base  one  upon  the  other  with 
their  bases  in  common  and  its  boundary  is  defined 
by: 


r  = 


sin(/r  /  4 )/  sin(3;r  /  4  -  6\  0  <  9  <  n  /  2,0  <  <f>  <  In, 
sin(;r/4)/sin(5/r/4  -  #),  nil  <  6  <  tz,0  <  <f>  <  2k, 


It  is  easy  to  see  that  f\  (o)  is  the  sum  of  an  infinite 
Fourier  series.  Figure  1  shows  the  reconstructions 
obtained  applying  the  procedure  described  in  Step 

k )-  h )  when  we  choose  Cl*  =  B,  c  =  1,  X  =  0. 1  or 
A  =  0.9,  n  =  1-X,  <J  =  1,  Lmax  =8,  e=5-10"3,  Ls  =8  or  Ls  =20, 
={4},E2  =  =  (sin(z;r/10),0,cos(/;r/10))r  e  dB,  i  =  1,2,...,1()| 

Z3  =  ^  =  (sin(/‘;r  /10)cos(j  7i  1 5),  sin(/;r  / 10)  sin  (7  k  1 5),cos(//z-/  io)f 

and  i  =  1,2,. ..,9,  y  =  l,2,...,10}u|(0,0,l)7’,(0A-l)7’| 

S4  =  (x  =  (sin((2/ -l)^-/20)cos((27  —  l)^r  / 1 0),  sin((2z  -  l);r/20) 
sin ((2  j  - 1)77 / 1 0), cos((2z  -\)n 1 20)r ,  i  =  1,2,. ..,10,  j  =  l,2,...,lo}  . 

We  note  that  the  reconstructions  shown  in  Figure  1 
are  very  similar  and  they  depend  essentially  only  on 

the  value  of  Ls  used  as  shown  by  the  reconstructions 

obtained  when  T  =  0.9  and  Ls  -  8  0r  Ls  -  20  t 
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X=0.9  1=8  >.=0.9  L  =20 

0  o 


Figure  1:  Furtive  smart  double  cone  and  three  reconstructions  versus  A  and  Ls 
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The  second  experiment  involves  a  smart  acoustically 
soft  (i.e.  Z  =  ®)  corrugated  sphere  whose  boundary 
is  given  by  r  =  /(6))  =  l-0.15sin2(4^),0<6,<^,0<^<2^- 
(see  Fig.  2))  that  pursues  the  goal  of  appearing  as  an 
acoustically  soft  (i.e  Z'a  =  ®)  double  cone  (see  Fig.  1 
up  left  comer)).  In  this  experiment  the  far  field  data 
used  for  the  reconstmction  are  synthetic  data 
obtained  solving  numerically  the  corresponding 
direct  scattering  problems  and  have  3-4  significant 
digits  exact.  We  have  chosen  the  initial  guess  of 
3Q*  assigned  in  Step  i\)  to  be  the  surface  of  the 
sphere  with  center  the  origin  and  radius  1.2, 
moreover  we  have  chosen  c  =  1,2  =  0.1  (Figure  2), 
2  =  0.9  (Figure  3),  //  =  i-2,^  =  i,zmax  =Lg  =8,e=5-i0“3, 
Ls  —  8,  Sj  =  {3},  S2  =  | at  =  (sin(/;r/10),0,cos(/;z710))r  edB,  i  = 

Z3  =  ^  =  (sin(z;r  /10)cos(y  7r  /  5),  sin(z';r  /  10)sin(y  7T  /  5),COs(i7T  /  io))r 

and  *  =  1,2,.. .,9,  j  =  l,2,...,10}u{(0,0,l)r,(0,0,-l)r|  and 

£4  =  {x  =  (sin ((2/  -  l);r/20)cos((2  j  -  l)/r/10),  sin ((2/  -  l)/r  /  20)  . 
sin ((2  j  - 1)71  / 1 0), cos((2/  -  \)7i /  20)r ,  i  =  1,2, ,...,10,  j  =  1,2, ...,lo}  . 
Figures  2  and  3  show  the  results  obtained  applying 
the  procedure  to  solve  the  inverse  problem  z0)-/7) 
when  2  =  0.1  (Figure  2)  and  2  =  0.9  (Figure  3).  As 
shown  in  the  figures  we  choose  different  values  of 
the  parameters  a  and  b  of  the  initial  guess  SQ  * . 
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a=-0.5  b=0.5  Q*=B]  2  a  =  -0.5  b=0.5  n*=B, 


Figure  2:  Masked  smart  obstacle  and  three  reconstructions  when  A  =  0.1 


a=-0.5  b=0.5  n=B,  2  a=-0.5  b=0.5  0*=B, 


Figure  3:  Masked  smart  obstacle  and  three  reconstructions  when  A  =  0.9 
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